Monitoring non-parametric profiles using adaptive EWMA control chart

To monitor the quality of a process in statistical process control (SPC), considering a functional relationship between a dependent variable and one or more independent variables (which is denoted as profile monitoring) is becoming an increasingly common approach. Most of the studies in the SPC literature considered parametric approaches in which the functional relationship has the same form in the in-control (IC) and out-of-control (OC) situations. Non-parametric profiles, which have a different functional relationship in the OC conditions are very common. This paper designs a novel control chart to monitor not only the regression parameters but also the variation of the profiles in Phase II applications using an adaptive approach. Adaptive control charts adjust the final statistic with regard to information of the previous samples. The proposed method considers the relative distance of the chart statistic to the control limits as a tendency index and provides some outcomes about the process condition. The results of Monte Carlo simulations show the superiority of the proposed monitoring scheme in comparison with the common non-parametric control charts.

www.nature.com/scientificreports/ follow the IC formulation. Note, the latter is not always a reasonable assumption because the OC model may not be easily determined (especially in complicated process, such as non-linear models) and it is often invalid in real-life practical applications. Moreover, it may increase the values of ARL 1 and non-compatible products (due to a high false alarm rate).
To overcome these problems, several non-parametric control charts have been extended for monitoring quality characteristics. In most of the non-parametric methods, a distribution-free statistic such as Mann-Whitney test is computed for the monitoring purpose 25 . Besides, several other approaches could be found in the related literature 26 . In spite of the fact that there are several non-parametric control charts for monitoring quality characteristics, few studies can be found for non-parametric profiles monitoring, especially in phase II applications. The following paragraphs provide a comprehensive review of the related literature in which the most focuses lie in the researches with different IC and OC models. Hereafter, similar to the main focus of this study, the profiles that have different IC and OC shapes are denoted as non-parametric profiles.
As the fundamental research, Williams et al. 27 extended five metrics in monitoring non-linear profiles in Phase I based on measuring the deviations from the IC situation. Although these methods created a parametric control chart, they can also be considered as non-parametric if the user does not assign a prespecified distribution of parameters. Considering this work, Zou et al. 28 developed a non-parametric multivariate EWMA (exponentially weighted moving average) statistic called NEWMA (non-parametric EWMA) for phase II profile monitoring. In this method, the coefficients and standard deviation parameters are firstly scaled and thereafter are taken into account as the previous samples to construct a MEWMA (multivariate EWMA) statistic. On the ground of simulations results of Zou et al. 28 's paper, it can be counted as the most fundamental control chart in phase II of non-parametric profile monitoring. Some other extensions of this fundamental paper could be found in Qiu and Zou 29 , Qiu et al. 30 , Hung et al. 31 , Chuang et al. 32 and Li et al. 33 .
Designing EWMA3 statistics proposed by Kim, Mahmoud and Woodall 7 considered as one of the most common approach in simple linear profiles (i.e. the EWMA3 statistic by Kim et al. 7 ), Zhang et al. 34 designed the corresponding non-parametric method. By consideration of changing the linear IC model to quadratic form, they compared the EWMA3 approaches with some hypothesis tests including LRT (likelihood ratio test), F and T 2 . The same approach has been utilized in Shang et al. 17 by considering the binomial and Poisson IC models. Similarly, Zhang et al. 34 and Zi et al. 35 developed a non-parametric control chart for monitoring linear profiles, where the Wilcoxon rank estimator was employed (instead of the common least squares approach) and the simulation results showed the favorable properties of this approach.
Because non-parametric methods usually deal with complicated profiles, wavelet transformation has been employed to simplify the complexity of monitoring and enhancing the identification ability of OC sources 36,37 . Also, to ease the computations of wavelets in non-linear profiles, PCA (principal component analysis) techniques has been suggested in Paynabar, Jin and Pacella 38 . Increasing the detection ability of control charts in nonparametric profiles with the aim of machine learning and ensemble learning could be found in Yeganeh et al. 39 in which more accurate control charts' techniques are used instead of conventional approaches. Jones et al. 6 provided some useful guidelines for the implementation of the non-parametric control charts in profile monitoring.
Most previous research on monitoring profiles focused on identical profile structure for OC and IC situations or parametric models; in other words, it is assumed that the OC models are restricted to have the same shape as the IC model. But this assumption may be violated in complicated real-life problems such as image and video processing, large-scale data, robotics, sensors' surveillance in a way that different shape between IC and OC profiles would be probable 6,17 .
Neglecting this issue leads to two main flaws in the monitoring procedure. First, the estimation of parameters does not fit with the usual methods like ordinary least square (OLS) and secondly, the next challenge occurs when the IC model is nonlinear. In the latter case, some specific estimators are required but they may not have a reasonable performance in case of changing the OC model shape. Zou et al. 28 discussed these situations and stated that definition of some non-parametric methods is necessary to avoid such similar problems. So, if there is an expectation of different IC and OC shape in a real problem, employing of non-parametric methods seems to be necessary to avoid false signals. It is also obvious that similar to parametric methods in phase II, quick shift detection ability for control charts (which is the main aim of this study) would be valuable.
Literature review of adaptive control charts in profile monitoring. In the previous literature of profile monitoring, there are several memory-type control charts especially for linear profiles. The idea of ranked set sampling was proposed by Riaz et al. 40 and Huwang et al. 41 in which the samples are taken as a memory to be ranked by some criteria. These rank-based monitoring schemes generates very sensitive statistic to small shifts. The Variable Sampling Interval (VSI) idea has also been studied in profile monitoring 14,[42][43][44][45] . The main idea of VSI is that the samples' interval is not constant and could be identified based on the OC occurrence probability.
As a general term, these approaches can be referred to as the adaptive control charts in which the sample size, statistic magnitude, interval between the sampling and other parameters are adjusted based on the estimation of the occurred shifts 14,40,42,[45][46][47] . To the best of authors knowledge, except Mohammadzadeh et al. 14 and Jeong et al. 47 , which have investigated adaptive methods based on the logistic and non-linear profiles, there are no other adaptive methods for complicated models such as the non-linear and non-parametric approaches. The main reason for this may be due to the weak performance and incompatibility with non-parametric limitations of adaptive approaches. Also, the complexity of parameter tunning could stir up some troubles in complicated profile structures.
The aims and innovations of the study. In the present study, a novel adaptive approach is developed for non-parametric profile monitoring in which there is currently no research based on the adaptive control www.nature.com/scientificreports/ charts in phase II SPC applications. The proposed approach, whose purpose is to improve the detection ability of NEWMA chart statistics, proposed by Zou et al. 28 , in term of ARL criterion, is developed based on a generally different scheme as compared to older adaptive mentioned methods. It considers the relative distance of NEWMA chart statistics to the OC situations as the major criteria in OC detection and then, converts the obtained relative distance to an EWMA statistics. To reach the best performance in phase II, i.e., minimum ARL 1 with a constant ARL 0 , a heuristic designing procedure of the parameters is implemented in this paper. The performance evaluation is conducted through extensive non-parametric simulations. In addition to the NEWMA approach as the base statistic, some other non-parametric and adaptive based competitors are used to show its superiority (when integrated with the heuristic design approach). Although the proposed method has been extended based on the NEWMA statistic, it is able to combine with other base statistics in phase II monitoring.
The main contributions of this study could be summarized as follows: • Proposing an adaptive control chart in profile monitoring.
• Reaching the optimum control chart parameters with a heuristic approach.
• Combination of the proposed adaptive approach with NEWMA statistic for monitoring non-parametric profiles. • Computing the detection ability of the proposed method in term of ARL in comparison with some nonparametric and adaptive conventional methods.
The remainder of this paper is organized as follows: the concept of NEWMA control chart is illustrated in "NEWMA statistic for non-parametric profiles". The proposed heuristic method of this study is discussed in "The proposed method". "Simulation results" provides the simulation results for performance evaluation and the corresponding sensitivity analyses is conducted in "Sensitivity analysis". A real-life example is used to illustrate the implementation of the proposed NEWMA chart in "Illustrative example". Finally, the concluding remarks and recommendations for the future research are presented in "Conclusion".

NEWMA statistic for non-parametric profiles
As we focus on the monitoring of profiles that could be well summarized by non-parametric regression model in phase II, the notations and formulations are extended based on the NEWMA scheme 28 as the most fundamental control chart in this field. Suppose that for the jth (j ≥ 1) random sample gathered over time, we have a set of observations entailing a n j -variate response vector (y j ) and p × n j matrix of explanatory (X j ) variables with sample size n j and for simplicity, it is supposed that the sample sizes are the same in all the profiles, so n j is replaced with n. Also the assumption of the same values in explanatory variable is very common in the related literature (i.e., using X instead of X j ) so we have X = (x 1 ,x 2 ,…,x n ). By these assumptions, the general IC model can be given by: where g 0 is the IC regression function and σ 0 2 is the IC parameter obtained from phase I or previous information about the process, which are assumed to be known in phase II. Considering G 0 = (g 0 (x 1 ), g 0 (x 2 ),…, g 0 (x n )), to monitor the above model at the jth sampling time with the vector of response variables Y j = (y 1j ,y 2j ,…,y nj ), Zou et al. 28 defined the following statistic: In Eq. (2), W is a symmetrical n × n smoothing matrix to be utilized in local linear estimator of response variables and obtained based on the explanatory variables, so it is fixed and constant in each new generated profile. To compute smoothing matrix, a symmetric probability density function and a bandwidth, denoted by h E are defined and each element of W are calculated based on them. For brevity, the formulas are not reported here but a MATLAB function for computation of W is provided by the authors and available upon request. For reaching optimum value of bandwidth, a relationship based on the explanatory variables was introduced in Eq. (10) of Zou et al. 28 . In that Equation, a constant parameter, denoted by c, is defined and the simulations are conducted based on different values of c.
For developing the MEWMA chart statistic, Zou et al. 28 first transformed the σ 2 j to a standard normal variable: In this transformation, φ −1 is the inverse of the standard normal CDF (cumulative distribution function) and φ can be approximated with the chi-square distribution in the way that: (1) www.nature.com/scientificreports/ Then, U j , denoted as (Z j ,σ j ) ′ , which is an (n + 1)-variate vector and the (n + 1)-dimensional symmetric covariance matrix is defined as = V 0 0 1 . The proposed EWMA charting statistic then could be calculated as: And finally, we have a positive chart statistic: where λ is the EWMA constant (here is equal to 0.2). Considering Eq. (6), an OC signal is triggered if Q j is located beyond the IC region Q j > L 2− , where L is assigned based on desired ARL 0 .

The proposed method
The major idea of the current study is to propose an adaptive approach that is established based on the tendency (relative distance) of obtained statistics to the OC situations. As mentioned in Aly et al. 48 and Haq and Khoo 49 , an adaptive control chart adjusts the chart statistic based on the shift size of the process in the current time. Following this idea, the proposed tendency or relative distance is computed according to the ratio of the statistic's distance from upper control limit (UCL) and lower control limit (LCL). So, the tendency of the jth sample (T j ) is defined as: In NEWMA approach, LCL is set at zero, so the above formula transforms to: Basically, the lower the charting statistics are, the greater tendency would be. For better understanding, Fig. 1 depicts the tendency values for three different NEWMA statistics when LCL and UCL are 0 and 2, respectively. For example, if the first NEWMA statistic is 0.5, T 1 is calculated as 3.
Naturally, such an OC situation in an NEWMA control chart will cause more samples with lower tendency values; in other words, the average of tendencies would be decreased in case an OC situation has been observed. So, the average of tendencies (T * j ) is firstly computed by the following equation and then, it is mapped to a specific rate for applying in an adaptive approach. www.nature.com/scientificreports/ The aim of this paper is to use the proposed tendency index in Eq. (8) as an auxiliary adaptive approach for adjustment of the final NEWMA statistic. To have a quicker signal detection in OC conditions in the proposed adaptive approach, the NEWMA statistic (Q j ) is increased (decreased) when the tendency index (T j ) is decreased (increased), respectively. The updating of the NEWMA statistic in the jth sample is done by using the tendency rate (defined by Tr j ) as follows: To have a more sensitive control chart against OC shifts, it is expected that Tr j to be lower than 1 when T* j has a large value (IC situation); note though, it can be greater than 1 when T* j has a small value (OC situation). Hence, a mapping function or relationship between Tr j and T* j will be required to use the proposed adaptive approach. This approach has been carried out in previous works by definition of a score function 48,50 , improved estimators 49 , warning limits 42,51 and so forth.
To compute Tr j , a completely linear function is not employed in this paper and a hinge function (called mapping function in this paper), as shown in Fig. 2, is utilized. A hinge function hinders very small rates causing very late signalling especially in large shifts so it is not allowed to have smaller tendency rates than lower limit (L L ) even the average of tendencies is greater than a predefined value b. In addition to a and L L , two parameters entailing B and upper limit U L have to be defined for obtaining the mapping function formulation. The dash line means that it is very unlikely to be a sample in that region.
Considering the known values for the a, b, L L and U L (designing of these parameters are discussed later), an OC signal will be triggered when Q* j is greater than UCL. This approach adapts the NEWMA statistic magnitude based on the previous samples in a way that the chart statistic becomes so large to easily reach an OC signal. It is shown through simulations that the proposed approach is able to gain tangibly lower ARL 1 values rather than NEWMA chart.
For better understanding, the procedure of ARL 1 simulations is depicted in Fig. 3. In this procedure, the OC profiles are generated to reach a signal in MaxIt times and the average of run lengths (signalling time) is reported as ARL 1 . Similar to ARL 1 , the standard deviation of run length (SdRL 1 ) is computed for each predefined shift in the process.
In this approach, the designing procedure involves selection a proper value for of a, b, L L and U L . Due to complex close form formulation, the control limits are adjusted based on the simulation in this paper; it is a common approach in previous profile monitoring studies 42,45,52 .
The procedure of designing comprises three main directions. First, the value of UCL is only related to the ARL 0 or equivalently, it is similar with the NEWMA approach. Considering the ration of sample in IC data generation for the coefficients i.e., a and b is the second one, and keeping the slope of mapping function at a predefined value is the third.
To this aim, 10,000 IC profiles are generated and T j (T 1 , T 2 ,…,T 9999 and T 10000 ) are computed. Then, the 45% and 95% quantile of these data are considered as proper values for a and b. These experimental values have been obtained after several investigations in different IC profiles. To adjust L L and U L , three following aims should be satisfied: • The final ARL 0 is supposed to be very close to its desired value. www.nature.com/scientificreports/ The rationale behind the first condition is that the proposed method tries to have greater statistics in OC situations. The second condition suggests a reasonable mapping function between Tr j and T * j . It has been obtained by several simulation studies in different IC profiles. Finally, the last one is the natural aim of each control chart in phase II application.
By three above conditions, we should be trying to reach the best values of L L and U L . In the above formula, L L is expected to be applied in IC condition while the chart statistic will become greater with U L in OC situation. In other words, the chart statistic is multiplied by a value which is less (greater) than one in large (small) tendencies (or IC (OC) situation), respectively. The second condition has been suggested by several investigations in different IC profiles and the last one is the main criteria for phase II applications. The following steps summarize the designing procedure.
• Step 1: Adjust UCL to reach desired ARL 0 . It is the same as the NEWMA approach.
• Step 2: Generate 10,000 IC profiles and store T 1 , T 2 ,…, T 9999 and T 10000 . • Step 3: Calculate the 45% and 95% quantile of these data as a and b.
• Step 4: Calculate L L and U L with regard to the three above conditions. • Step 5: Obtain the mapping function formula.

Simulation results
To compare our proposed method, three different simulation setups in which the IC model respectively has the polynomial, exponential and linear form are adjusted here, based on the non-parametric approaches simulated in Zou et al. 28 and Zhang et al. 34 . The IC model of these scenarios are defined in Table 1 for the jth generated profile (i = 1, 2,…, n).
Note that in all the models σ 0 = 1. The competitive methods for the first and second scenarios are NEWMA, NM (naïve MEWMA) and PM (parametric MEWMA). All of these methods have been defined in Zou et al. 28  (2) and for brevity, we only consider the last situation (c = 2) in the first and second scenarios. To start the designing procedure (discussed in previous section), the values of UCL are considered as the values reported in Tables 2 and 4 of Zou et al. 28 , i.e. it is computed with L 2− . In scenario 3, the competitive methods are KMW (it is the same as EWMA3 approach in Kim et al. 7 ) and NEWMA with c = 1 (it was not simulated in Zhang et al. 34 and the results have been obtained based on our simulations). To reach ARL 0 equal to 200 in NEWMA approach, L is set to 19.58 and 20.75 for n = 13 and 25 (the initial values of UCL in our procedure). It is worthwhile to mention that Zhang et al. 34 did not propose a clear manner for determination of control limits of KMW. But we speculate that they adjusted each separate chart to reach the same ARL 0 . Note that simulations were not conducted for the case of n = 5, this is due to the inability of NEWMA method to signal in some shifts when dealing with small sample sizes (for more details, see Remark 1 in Zou et al. 28 ).
We tried to define similar models as those in Zou et al. 28 and Zhang et al. 34 ; however, due to some inconsistency in the results of second scenario in Zou et al. 28 , there are some differences in this case. The OC models in scenario 1 are defined as: In this scenario, six different OC shifts have been employed based on Table 2.
In the second scenario, the OC models are: In this scenario, six different OC shifts have been employed based on Table 3.
In the third scenario, Zhang et al. 34 supposed that the IC model changed to quadratic form by adding the γ (x * 2 i − η) term to the IC model in the way that γ = δσ. Note that the EWMA constant is assumed equal to 0.2 in all scenarios.
Considering ARL 0 equal to 200, the values of UCL in each scenario based on the Zou et al. 28 's method (NEMWA approach) are reported in Table 4. Note that the results of first and second scenario are the same as Tables 2 and 4 in Zou et al. 28 . Also, the design parameters for each scenario were gathered in Table 4.
In the following subsections, our proposed method is compared with other competitors based on the ARL 1 criteria in each of the predefined scenarios.     Table 5. These results reveal that the proposed method (ANEWMA) is effective for all of the shifts. For large shifts in the process parameter, its performance is roughly better than competing methods. As can be seen in the results of model (II) and (III), the more complicated the OC model is, the better the outperformance.    Table 3 (it is the same as Table 3 in Zou et al. 28 ). We may be asked about the changing of the third OC model from Zou et al. 28 's proposed model y ij = 1 1+β 1 x β 2 i + ε ij . This is due to our inability to produce the same results with Zou et al. 28 's paper. We even re-simulated the results of the NM method (PM was not able to detect these types of shifts) to make sure there was nothing wrong with our simulations, but unfortunately those results were also not obtained. It may be possible that the third model was written incorrectly in that paper.

IC model
The results of ARL 1 with respect to the above OC models are gathered in Table 6. These results reveal that the proposed method is effective for all of the shifts and the same conclusions with scenario 1 can be drawn and only in few cases, ANEWMA is not the best approach.
The simulated ARL 1 values for scenario 3. The results of ARL 1 for these types of shifts for n = 13 and 25 are reported in Tables 7 and 8, respectively. In these tables, the results of EWMA3 approach are shown with KMW notation.
The superiority of the NEWMA to KMW is predictable due to the linear structure of the KMW control chart. Also, NEWMA has better performance than the other competitors reported in Zhang et al. 34 when n = 13. However, when n = 25, the ANEWMA outperformed the other competitors in small shifts only, while the NEWMA is the best for the detection of moderate and large shifts. Since the ANEWMA naturally improves the detection ability of NEWMA, it gives logically the best performance for the detection of all the shifts.

The simulated SdRL 1 values for scenario 3.
To evaluate control charts in phase II, SdRL 1 could be used in addition to ARL 1 53-55 and such a proper approach is expected to have lower SdRL 1 values. Figure 4 illustrates the SdRL 1 values for ANEWMA, NEWMA and KMW schemes. We can see that the proposed approach was able to reduce SdRL 1 for ANEWMA chart in comparison with NEWMA chart. The same results have been generated by other scenarios but they have been omitted for the sake of brevity.

Sensitivity analysis
In this section, two different simulations are provided to show the superiority of the proposed method. In the first simulation, two conventional adaptive approaches are prepared for monitoring non-parametric profiles. In the first approach, following the Zou et al. 23 and Mohammadzadeh et al. 14 , VSI scheme is implemented and the proposed adaptive method in Haq 46 is the second one. These methods hereafter called VSINEWMA and HQNEWMA are incorporated with NEWMA statistic to have a fair comparision. For brevity, the details of designing have not been provided here.
To show the effectiveness of the proposed adaptive method, the NM statistic 28 are combined with the proposed adaptive approach [denoted with ANM (Adaptive NM)] in the second simulation, the purpose of which is to compare the detection ability of ANM and NM. To this aim, tendency rate based on Eqs. (9) and (10) is multiplied   www.nature.com/scientificreports/ by NM statistic instead of NEWMA statistic. The designing procedure is the same as the ANEWMA and for the brevity the details are neglected. Also, the results of the first scenario are only reported due to page limit.
Comparison of ANEWMA with other adaptive approaches. Figure 5 reports the ARL 1 values of the first scenario for model (I) and (II) when n is 20 and 40. Note that the ANEWMA results have been gathered from Table 5. The superiority of the proposed method over two adaptive competitors is obvious from the following simulations. Except the shift type (v) in model (I) when n = 20, ANEWMA has the minimum ARL 1 value in all of the simulations. Also, VSINEWMA is nearly better than the HQNEWMA in these simulations.
The effect of proposed adaptive approach in combination with NM statistic. Figure 6 reports the ARL 1 values of the first scenario for model (I) and (II) when n is 20 and 40. In all of the shifts, ANM had the better performance over NM control chart which is an indicator of the capability of the proposed adaptive approach. Also, ANEWMA outperformed ANM for all the shifts which is due to the better detection ability of NEWMA, compared to NM. The supririty of the proposed method over two adaptive competitors is obvious from the following simulations. Except the shift type (v) in model (I) when n = 20, ANEWMA has the minimum ARL 1 value in all of the simulations. Also, VSINEWMA is better than HQNEWMA in these simulations. The above two simulations revealed two main findings. First, the proposed method not only outperformed non-parametric and parametric control charts but also had the superiority over conventional adaptive methods. As the second finding, it may be possible to reach better performance than the Haq 46 's method by using the optimisation of the parameters regarding the non-parametric profiles; however, it is not in the scope of this paper and would be valuable for the future research.

Illustrative example
A real case is provided in this section to show the practical application of the proposed scheme. It was selected from the semiconductor device fabrication in which there is a deep reactive ion-etching process. The aim of this process is to check an etched wafer by an electron microscopy scanner in scientific labs. Employing of this wafer is needed to run complicated mechanical-chemical reactions for a complex automotive system 23,24,28,56 . In these reactions, one key factor is to determine the profile of a trench in the system because it can change the outcome of the downstream operations and also the final quality of the product.
To establish a suitable monitoring approach for such a trench profile, Zou et al. 23 suggested a new method based on the profile monitoring control charts. By using the Scanning Electronic Microscope (SEM) data rather direct measurements, they were able to reach at a proper monitoring scheme.
They suggested an IC polynomial model (Eq. 11) for this problem with the following transformation z i = x i y ij = 1.55 + 0z i + 0.62z 2 i + ε ij , i = 1, 2, . . . , 11; j = 1, 2, . . . , x i = −2.5 : 0.5 : 2.5, ε ij ∼ N(0, 0.16), www.nature.com/scientificreports/ where the dependent and independent variables are the shape of the profile. Due to complexity of getting real data, the new samples were generated by simulations in Zou et al. 23 . To have an ARL 0 equal to 370, UCL NEWMA was set at 2.01 (L = 18.00) when c = 1.5 28 . So, the design parameters (as described in "The proposed method") are a = 2.3, b = 4.6, L L = 0.46, and U L = 2.85. The mapping function's intercept and slope were obtained as 5.24 and − 1.04 and 1.41 by the proposed designing scheme. By changing the standard deviation to 0.48 from 0.4 (it is equal to a shift magnitude of 1.2σ), the OC profiles were generated. Table 9 gathers the details of the 15 OC generated profiles. As can be seen, NEWMA was able to trigger an OC signal in the 15th sample while ANEMWA only needed 7 samples to detect this shift.
For better understanding, the computations of Tr j (in the second and last sample) are illustrated here.  Fig. 7 with green and black lines, respectively. The values of Tr j are shown in red.

Conclusions
In this paper, an effective adaptive approach for improvement of the non-parametric profile monitoring in Phase II was proposed. The adaptive approach was combined with the NEWMA control chart proposed by Zou et al. 28 , to enhance the OC detection in Phase II monitoring. To compute the control limits, a simulation-based designing procedure was proposed based on the Monte Carlo theory. The results showed that the proposed method not  www.nature.com/scientificreports/ only improved the performance of NEWMA approach in phase II, but also outperformed some of the existing control charts, such as NM, EWMA3 and other adaptive methods, for most shifts in terms of ARL criterion. In these simulations, three different scenarios (IC models) entailing polynomial, exponential and linear models have been investigated to show the robustness of the proposed scheme. In the case study, an example from the semiconductor device fabrication was implemented in which remarkable ability of the proposed method for the detection of OC shifts was demonstrated. For future research, we suggest the development of a control chart based on the combination of this approach and other adaptive control schemes. Also, obtaining the exact distribution of the chart statistics and control limits for adaptive control charts would be a novel idea.